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We provide a quantitative description of the structure of edge states in split-gate quantum wires in 
the integer quantum Hall regime. We develop an effective numerical approach based on the Green's 
function technique for the self-consistent solution of Schrodinger equation where electron- and spin 
interactions are included within the density functional theory in the local spin density approximation. 
The major advantage of this technique is that it can be directly incorporated into magnetotransport 
calculations, because it provides the self-consistent eigenstates and wave vectors at a given energy, 
not at a given wavevector (as conventional methods do) . We use the developed method to calculate 
the subband structure and propagating states in the quantum wires in perpendicular magnetic field 
starting with a geometrical layout of the wire. We discuss how the spin-resolved subband structure, 
the current densities, the confining potentials, as well as the spin polarization of the electron and 
current densities evolve when an applied magnetic field varies. We demonstrate that the exchange 
and correlation interactions dramatically affect the magnetosubbands in quantum wires bringing 
qualitatively new features in comparison to a widely used model of spinless electrons in Hartree 
approximation. 
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I. INRODUCTION 

Transport properties of quantum dots, antidots and 
related structures are affected by the nature of current- 
carrying states in the leads connecting these structures 
to electron reservoirs. In sufhciently high magnetic fields 
the current-carrying states arc the edge states propagat- 
ing in a close vicinity to the sample boundariesi. A de- 
tailed information on the structure of the edge states rep- 
resent a key to the understanding of various features of 
the magnetotransport in the quantum Hall regime. 

A quantitative description of the edge states for 
the case of the gate-induced confinement of the high- 
mobility two-dimensional electron gas (2DEG) was given 
by Chklovskii et al^, who provided an analytical solu- 
tion for the positions and widths of the compressible and 
incompressible strips arising in the 2DEG due to the elec- 
trostatic screening. In the compressible regions, the Lan- 
dau bands are pinned at the Fermi energy Ep. This leads 
to a metallic behavior when the electron density is redis- 
tributed (compressed) to keep the electrostatic potential 
constant. In the incompressible regions, where the Fermi 
energy lies in the Landau gaps, all the levels below Ep 
are completely filled and hence the electron density is 
constant (which is consistent with the behavior of the 
incompressible liquid) . 

A number of studies addressing the problem of 
electron-electron interaction in quantum wires be- 
yond the electrostatic treatment of the edge states of 
Chklovskii et al^ have been reported during the re- 
cent decado'^i^i^'^i'^i^^i^°i"i^^i"i^^. The many-body as- 
pects of the problem have been included within Thomas- 
Fermi^, Hartree-Fock^i^, screened Hartree-Fock^, and 
the density functional theorjt^. The full quantum- 
mechanical calculations based on the self-consistent solu- 



tion of the Schrodinger equation have been done within 
the HartreoifliiiiiSii^ and the density functional theorjiii 
approximations. 

A particular attention has been devoted to investiga- 
tion of the spin polarization effects in edge states^i^iSiL&. 
For example, Dempsey et al^ have shown that for a suf- 
ficiently smooth confining potential, spin degeneracy of 
the outermost edge state is lifted and two spin channels 
become spatially separated. The interest to the spin- 
related effects in quantum wires is also motivated by 
significant current activity in semiconductor spintron- 
ics, that utilizes the spin degree of freedom of an elec- 
tron to add the additional functionality to electronic de- 
vices. A number of proposed and investigates devices 
for spintronics applications operates in the edge state 
regimaiSi^, which obviously requires a detailed knowl- 
edge of the spatial dependence of the spin-resolved states 
in the quantum wires. Edge states have also been pro- 
posed as one-way channels for transporting quantum 
informatiornti. The knowledge of the spin/charge struc- 
ture of the current carrying states is also essential for 
numerical simulation and modelling of the magnetotrans- 
port in quantum dots and related structures (Note that 
such the modelling is often done utilizing single-electron 
wave functions in the leads disregarding the spin/many 
electron effect siSiiSiSS) . In order to obtain such the in- 
formation on quantum-mechanical propagating states in 
quantum wires, one has to solve the Schodinger equa- 
tion incorporating the exchange interaction to account 
for the spin effects. In should be noted that the stud- 
ies reported so far are often limited to some strictly 
integer filling factors^, or utilize Thomas-Fermi-type 
approaches^ or perturbative techniquo^, where the re- 
quired information concerning the quantum-mechanical 
wave functions is not available. Moreover, the quantum- 
mechanical effects associated with the finite extension of 
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the wave function (not included in e.g. Thomas-Fcrmi 
approach) can play a decisive role for the quantitative 
description of the edge states. For example, Suzuki 
and Andoi^ have demonstrated (in a model of spinless 
electron) that the predictions of Chklovskii et al. and 
Thomas-Fermi models regarding the existence and the 
size of the compressible/incompressible strips are in qual- 
itative disagreement with the self-consistent modelling 
based on the Schrodinger equation in the regime when 
the estimated width of the strips is smaller than the ex- 
tend of the wave functions. 

The purpose of the present article is two- fold. First, 
we perform a detailed self-consistent solution of the 
Schrodinger equation incorporating spin/many-body ef- 
fects in quantum wires. We discuss how the spin-resolved 
subband structure, the current densities, the confining 
potentials, as well as the spin polarization of the elec- 
tron and current densities evolve when an applied mag- 
netic field varies. We demonstrate that the exchange and 
correlation interactions dramatically affects the magne- 
tosubbands in quantum wires bringing qualitatively new 
features in comparison to a widely used model of spin- 
less electrons in Hartrce approximation. In the present 
study we limit ourself to the regime when more than one 
spin-resolved state can propagate in the wire, i.e. the fill- 
ing factor V > 1 (The filling factor v = n/nB = 27rZ^n, 
where n is the sheet electron density, tib = eB/h is the 
number of states in each Landau level per unit area, and 

= •\/^ is the magnetic length). 

Second, we present a detailed description of the devel- 
oped method based on the Green's function technique for 
the calculation of the subband structure and propagating 
states in the quantum wires in the magnetic field. This 
method is numerically stable, and its efficiency is related 
to the fact that calculations of the wave functions and 
wave vectors are reduced to the solution of the eigenvalue 
problem (as opposed to the conventional methods that re- 
quire less efficient procedure of the root searchingiSiiiii^) . 
The major advantage of the present method is that it 
can be directly incorporated into magnetotransport cal- 
culations, because it provides the eigenstates and wave 
vectors at the given energy, not at a given wavevector 
(as the conventional methods do). Besides, the present 
method calculates the Green's function of the wire, which 
can be subsequently used in the recursive Green's func- 
tion techniqueiSiSi widely utilized for magnetotransport 
calculations in lateral structures. 

In order to incorporate the spin/many-body ef- 
fects into the Scrodinger equation we use the den- 
sity functional theory (DFT) in the local spin-density 
approximation^^. The choice of the DFT is motivated, 
on one hand, by its efficiency and simplicity in the practi- 
cal implementation within usual self-consistent formula- 
tion introduced by Kohn and ShamS^, and, on the other 
hand, by its success in the reproduction of the electronic 
and spin properties of the low-dimensional structures in 
comparison to the exact diagonalization and quantum 
Monte-Carlo calculations, as well as experiments (for a 



review, see24). For example, Ferconi and VignaleS^ find 
that the accuracy of the DFT for the energy and density 
of few-electron quantum dots yields the accuracy better 
than 3% in comparison to the exact results. An excel- 
lent agreement between DFT and the variational Monte- 
Carlo results for the chemical potential and the addition 
spectra of the rectangular quantum dot was reported by 
Rasanen et alm^. 

Within the local spin density approximation the ex- 
change and correlation potentials are calculated using 
a parameterization of the functional for the exchange 
and correlation energy Cxc- The latter is usually obtained 
on the basis of quantum Monte Carlo calculationaSLS^ 
for corresponding infinite homogeneous system. In the 
present paper we use the parameterization of Tanatar 
and Cerperly (TC)^. This parameterization is valid for 
magnetic field when i/ > 1 , which defines the range of ap- 
plicability of our results. (Various parameterizations for 
e^c for strong fields v < 1 as well as different interpola- 
tion schemes between the low and the strong fields are re- 
viewed in Refsi2^i2S). Note that the DFT was used for the 
description of the spin polarization of the edge states in 
quantum wires in the integer Hall regime within Thomas- 
Fermi approximation^, as well as for the treatment of 
spinless edge states in the Kohn-Sham scheme based on 
the solution of the Schrodinger equationi^. The density 
functional theory within the Thomas-Fermy approach 
was also applied for the description of the edge channels 
in the quantum wire in the fractional Hall regime, where 
the parameterization of exc incorporated the additional 
gaps that open up at the fractional filling factors^. 

The paper is organized as follows. In Section II we 
present a formulation of the problem, where we define 
the geometry of the system at hand and outline the self- 
consistent Kohn-Sham scheme within the LSDA approx- 
imation. In Section III we provide a detailed description 
of or method based on the Green's function technique, 
and Section IV presents the major results and their dis- 
cussion. The conclusions are given in Section V, and 
Appendix presents some technical details of the calcula- 
tions. 



II. FORMULATION OF THE PROBLEM 

We consider an infinitely long split-gate quantum wire 
in a perpendicular magnetic field. A schematic layout 
of the device is illustrated in Fig. ^(a')- The distance 
between gates is a, the distance from the surface to the 
electron gas is b (we disregard the spatial extension of 
the electron wave function in the ^-direction) . The donor 
layer with the donor density pd has the width d and is 
situated at the distance c from the surface. The external 
electrostatic confinement potential can be written in the 
form 

chottky 1 (1) 
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FIG. 1: (Color online), (a) A schematic layout of a split-gate 
quantum wire in a perpendicular magnetic field, (b) Solid 
line: the calculated electrostatic potential Vconf{y) for the 
quantum wire with a = 500 nm, & = 60 nm, c = 14 nm, d = 36 
nm, pd = 6 • 10^^ Hi"'', and |Vb| = 0.2 V. A dashed line shows 
the corresponding parabolic confinement VLj(j/) = m* (ujy)'^ /2 
(with huj = 2meV) often used to approximate the external 
electrostatic confinement. 



a perpendicular magnetic field, = Bz, 

= Ho + Vconf (y) + K// (y) + 9^^bBc7, (4) 

where Hq is the kinetic energy in the Landau gauge, A = 
(-By, 0,0), 
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The last term in Eq. Q accounts for Zeeman energy 
where Hb = eh/2me is the Bohr magneton, a = ±i de- 
scribes spin- up and spin-down states, t , ij and the bulk 
g factor of GaAs is g ~ —0.44. 

The effective potential, Vfj:j{y) within the framework 
of the Kohn-Sham density functional theory reads^SiSiSi, 



V^ffir) = VHiy) + V,^M. 



(6) 



Vh (r) is the Hartree potential due to the electron density 
n{y) ~ n'^iy) (including the mirror charges). 
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The exchange and correlation potential Vxdy) — Vx{y) + 
Vc{y) in the local spin density approximation is given by 



where Vg (y) and Vd are respectively potentials due to the 
and the donor layer—, and Vscho is the Schottky 
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barrier. 
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(3) 



with Vq being the (negative) applied gate voltage, and 
Er being the dielectric constant. The Schottky potential 
is chosen to be Vschottky = 0.8 eV, which is appropriate 
for GaAs. The external electrostatic confinement poten- 
tial is shown in Fig. ^ (b) for a representative quantum 
wire with parameters typical for an experiment. Figure^ 
(b) also illustrates the corresponding parabolic potential 
V{y) = ^ (wy)^ often used to approximate the electro- 
static confinement in the split-gate wires, where m* is 
the effective electron mass (m* = 0.067me for GaAs). 
The wire is described by the effective Hamiltonian in 



=^Ke,e(n,C(y))} 



(8) 



where ({y) = is the local spin-polarization. As we 

mentioned in Introduction, in the present paper we use 
the parameterization of Tanatar and Ccrpcrly (TC)Si; 
for the sake of completeness, the explicit expressions for 
^^y) and Vc{y) are given in Appendix (see also Ref. 



III. CALCULATION OF THE ELECTRON 
DENSITY AND EDGE STATES IN QUANTUM 
WIRES. 

In order to calculate the self-consistent electron den- 
sities, wave functions and wave vectors of the magneto- 
edge states as well as corresponding currents, we use the 
Green's function technique. A detailed account of the 
major steps of the calculations is presented in this sec- 
tion. 
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A. Hamiltonian in the mixed energy-space 
representation. 

Numerical computation of the self-consistent electron 
densities and other quantities of interest requires the dis- 
cretization of the Hamiltonian (0J. Introduce a numerical 
grid (lattice) with the discrete variables m, n according 
to x,y ^ ma,na, where a is the lattice constant. The 
computational domain consists of Ns sites in the trans- 
verse n-direction (the wire is infinite in the longitudinal 
m-direction). Discretization of the continuous Hamilto- 
nian Q gives a standard tight-binding Hamiltonian with 
the magnetic field included in the form of the Peierls 
substitutionSi, 




n)} a. 



(9) 



where 



t {'^m,n'^m,n+l + e *''"a^ „a,„+i ^„ -|- h.c. } >, 



Vin) = Vconfin) + V^Jn) + gii^Ba (10) 



is the total confining potential, the hopping element 
t = /2m*a^ ^ the site energy eg = 4i, and q = eBa^/h; 

„ and am,n denote the creation and annihilation op- 
erators at the site (m, n) . The translational invariance in 
the longitudinal direction dictates the Bloch form for the 
propagating states in the quantum wire, 



(11) 



n=l 



where the index a corresponds to the a-th Bloch state 
with the wave vector k'^ and the transverse wave func- 
tion tlj^{n). In Eq. 1)11(1 the wave function tpai''^) cor- 
responds to the real space representation. To facili- 
tate the numerical calculation, it is convenient to ex- 
pand the wavefunctions over the transverse eigenstates 
(modes) of a homogeneous wire of the width of Ng sites, 

^^.(„) ^ . rjZo^r. JLUL. 



N, + l 
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N 



(12) 



in the real space representation and N in the "energy" 
representation), passing to the "energy" representation 
greatly enhances the computational speed. For example, 
for the wire of the width of 0.5 /xm, it is sufficient to 
use N Ki 50 modes to achieve a good convergence of the 
results with respect to the mode number. At the same 
time, in the real space representation, Ns = 500 (with the 
lattice constant a = 1 nm), which makes computations 
rather impractical. 

Passing from the real space representation to the "en- 
ergy" representation in the transverse direction we arrive 
to the Hamiltonian in the mixed energy-space represen- 
tation (i.e., in the real space representation in the longi- 
tudinal m-direction and "energy" representation in the 
transverse n-direction)i&, 

{N N 



^ L + R + 1 



(13) 



where 



= 2t 



2t cos j^j-^ are the eigenvalues of the 
transverse motion corresponding to the eigenfunctions 
(j)j{n); the creation and annihilation operators in the 
mixed space-energy representation are related to the real 
space creation and annihilation operators according to 



«mj = Y.nLi't>3{'n)a'tn,n^ a,nj ~ El^i^iWa™:"- The 
matrix elements of total confining potential and the hop- 
ping matrix elements read are given by 



tfr ^tY^Hn)e^'^-4>r{n), tf^, = (t« )*■ 



Note that Hamiltonian (|13|l has nearest-neighbor cou- 
plings in the longitudinal m-direction (described by two 
last terms in Eq. H13|) ). In the transverse ("energy") di- 
rection, the magnetic field couples all states j on slice m 
to all states jf on neighboring slices m -I- 1 and m — 1. 
The Bloch wavefunctions (|11|) in the mixed space-energy 
representation read 



N 



|0) 



(15) 



where the expansion coefficients "0^ j can be interpreted 
as the wavcfunction in the "energy" representation in the 
space of the transverse eigenstatesi^. Note that Eq. H12() 
corresponds to a conventional sin-transformation, whose 
inverse transform is given by the same equation^. The 
summation in Eq. (|12|l runs over 1 < j < N , with 
TV = Ns- In practice, however, it is sufficient to limit 
the summation to much smaller number of modes, with 
TV <C iVs ■ Because the speed of the method is determined 
by the dimension of the matrices (that is given by Ns 



B. Bloch states of a quantum wire in magnetic field. 

Define a retarded Green's function of the Hamiltonian 
iJ in a standard wa\i2iiM, 



iE-H + ie)G= 1, 



(16) 



where 1 is a unitary operator. Calculate first the Green's 
function g"^ corresponding to a single slice (see Fig. |21 
(a)). The Hamiltonian of m-th single slice reads 
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(a) Green's function of a single slice 
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(b) Semi-infinite wire 
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(c) Green's function of an infinite wire 
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FIG. 2: (Color online). Graphical illustration of the calcula- 
tions of the Green's function for a single slice (a), and for an 
infinite wire (b), (c). 
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(17) 



(Note that a single slice is not coupled to its neigh- 
bors, and hence two last terms in Eq. (|13|l are absent 
in Eq. (|17|l ). Using this operator in the definition of 
Green's function Hlt)|) . and taking the matrix elements 



(0|a, 



we arrive to the NxN system of linear 



equations for the matrix elements of the Green's function 

,4- 



of a single slice gj„j, = 



N 



(18) 



j"=i 



Note that because of the translational invariance, we have 
dropped index m in the definition of the matrix element 
of the Green's function of a single slice. 

Knowledge of the Green's function of a single slice g"^ 
allows one to find the Bloch states of an infinite wire. The 
eigenvectors {ipja} a-nd eigenvalues {fc^} are determined 
by the eigenequatior 



as 



Ak 



■0" 



if 



(19) 

where the matrixes t^,t^ and g'^ have matrix elements 
given by Eqs. and ifT^ respectively, and tp'^ is the 
column vector composed oi ip'j ^ I < j < N (Eq. ^^). 
Here and hereafter we use Greek indexes a, a' for Bloch 
states of the wire, and Roman indexes j, j' for the ba- 
sis set of the transverse eigenfunctions {4>j{n)}. Equa- 
tion H19|l has 2N eigenvalues k'^, 1 < a < N , which can 



be real or complex, describing respectively propagating 
and evanescent states (Here and hereafter k'^ is given 
in units of a^^ and the group velocity v is in units of 
a). The eigenvalues corresponding to right propagating 
states (w = dE/dk > 0, 5(A:) = 0) and states decay- 
ing to the right (5(fc) > 0) we denote by k^ with cor- 
responding eigenstates VJa ■ Correspondingly, k^ and 

■0Jq, stand for left propagating states {v = dE/dk < 0, 
3(/c) = 0) and states decaying to the left (3(fc) < 0). 
Sorting right- and left-propagating eigenstates can be 
easily done by calculating their group velocityi^ 



N 



(20) 



Passing to the real space representation for the wave 
functions in the above expressions and using the 
quantum-mechanical particle current density for a-th 
Bloch state calculated in a standard way for a tight- 
binding latticei^, 



If 

J- (n,i?) = - sin (qn + fc^) 1^-2(^)12, 



(21) 



the group velocity H20() can be expressed as the total par- 
ticle current of a-th Bloch state. 



(22) 



To conclude this section we note that a direct calculation 
of the eigenvectors fc^ and eigenfunctions ^'J/^ by substi- 
tution of Eq. (|15|1 into Schrodinger equation and calcula- 
tion of the roots of the corresponding determinant is pos- 
sible (see e.g., Refs. ITollTlllTa) . However, the procedure 
used here is more efficient as the solution of the eigen- 
problcm H19|l is numerically faster and less demanding 
than the root searching. Besides, an important advan- 
tage of the present method is that it can be directly in- 
corporated into magnetotransport calculations, because 
in contrast to the root searching method, the present 
technique provides the eigenstates and wave vectors at 
the given energy, not at a given wavevector. 



C. Calculation of the local electron density. 

The diagonal elements of the total Green's function of 
an infinite wire in the real space representation give the 
local density of states (LDOS) at the site r = (to, n)^. 



p'^{v,E) = --^[G%v,v,E)] 



(23) 



The LDOS p'^(r, r, E) can be used to calculate the local 
electron density at the site r, 

/>OC> 

n%v)^ dEp^{v,E)f{E-EF), (24) 

JVh 
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FIG. 3: A typical integration contour used in the calculation 
of integral II24I I. Dots indicate the poles of the Fermi-Dirac 
distribution function in the upper complex plane at J? [E] = 
Ef, Q [E] = (2m + l)7rfcT, m = 0, 1, 2, . . . . 

where f{E — Ep) is the Fermi-Dirac distribution func- 
tion and the lower limit of integration Vt corresponds to 
the bottom of the total confining potential. Note that 
p^iv, r, E) is a rapidly varying function of energy diverg- 
ing as ^ (iJ — Ea) ' when E approaches the threshold 
subband energies Ea. Because of this, a direct integra- 
tion along the real axis is rather ineffective as its numer- 
ical accuracy is not sufficient to achieve convergence of 
the self-consistent calculation of the electron density . We 
therefore calculate integral (|24|l by transforming the inte- 
gration contour into the complex energy plane 3 [E] > 
where the Green's function is much more smoother than 
on the real axis. (Note that all poles of the Green's func- 
tion H23|l are in the lower half-plain 3 < 0). A typi- 
cal contour used in the integration avoiding poles of the 
Fermi-Dirac function f{E — Ep) is shown in Fig. O Wc 
calculate the diagonal elements of the total Green's func- 
tion G"^(r,r, i?) as follows. We start from a semi-infinite 
quantum wire and calculate its surface Green's function 
r (i.e. Green's function for the boundary slice m = 1), 
sec Fig. 111(b). The right and left surface Green's func- 
tions r*" and r' (corresponding to a semi-infinite wires 
open respectively to the right and left) can be written in 
a matrix formiS*^ 

Tit"' = -^tKt-^t" (25) 

where the matrix elements {^t^~'')ja = 
{Kt^~^)aa' = exp ^fc^^' Saa'- Wc then connect this 
semi-infinite wire to the second semi-infinite wire to form 
an infinitely long quantum wire as shown in Fig. [21 (c). 
The total Green's function can be calculated with the 
help of Dyson equationSLM 

G = G^ + G°UG, (26) 

where correspond to the "unperturbed" structures 
(the left and right semi- infinite wires), and the operator 
U describes the interaction between them, 

N 

^ = - E + ^fj'«u«o.,'] , (27) 



(see Eq. (|13|) '). Using Eq. (|26() to calculate the matrix 
element {G'^)jji ~ (0|aij G"^ Oj'^^, |0), we obtain Green's 
function for the slice m = 1 (Note that because of the 
translational invariance in the m-direction, the calculated 
Green's function is the same for all slices), 

G^ = (i-r;t^r;,t^)"'r;, (28) 

where 1 is the unit matrix. Note that Eq. H28() 
gives Green's functions in the "energy" representation 
of the space of the transverse eigenstates. To obtain 
the Green's function in the real space representation 
needed to compute the electron density n'^(r) we 
perform a standard change of the basis, G'^(n,n, £') = 

D. Self- consistent calculations 

Iteration procedure. We calculate magneto-edge states 
and electron densities in a quantum wire in a self- 
consistent way, when on each iteration step a small part 
of a new potential H1(J|I is mixed with the old one (from 
the previous iteration step), 

Vr+i{n) = (1 - e)Vr{n) +eV^,{n), (29) 

e being a small constant, ~ 0.1 — 0.01. Using this input 
potential we, for a given energy E, solve the eigenproblcm 
(|19|l to find the Bloch states in the quantum wire. (Note 
that energy E is chosen in the complex plane as shown in 
Fig. (PJ). We then use the obtained results to calculate 
the Green's function G'^(n, n, E) according to Eqs. 
()28|l . The integration of the Green's function (|24(l gives 
the electron densities n" (r) , which are subsequently used 
to compute the new total confining potential H10|) . It is 
typically needed ~ 1000 iteration steps to achieve our 
convergence criterium W-^i) —n\^\/ {v}-^j^ +n\^) < 10^^, 
where n\jj is the one dimensional electron density ni/j = 
J n(r) dy on i-th iteration step. 

Adjustment of the Fermi energy. When the same fixed 
Fermi energy Ep is used for different magnetic fields -B, 
the calculated self-consistent one dimensional electron 
density changes as B varies. Depending on a particular 
realization of a quantum wire, one might need to adjust 
Ep for each B in order to keep the total electron den- 
sity fixed, as magnetic field does not change the electron 
density in the system. However, in a typical experimen- 
tal situation when a long quantum wire is connected to 
a 2DEG™, the Fermi energy in the reservoirs (not the 
electron density in the wire) is fixed. Because of this, in 
all calculations reported in the paper we keep Ep fixed 
(we set Ep = 0). 

Note that we have also performed calculations where 
Ep was adjusted to keep the electron density con- 
stant. All the results obtained in this case (in particu- 
larly, the density and current spin polarizations) are qual- 
itatively and quantitatively similar to those obtained in 
the case when Ep is adjusted. 
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Block states, subband structure and current density. 
Having calculated the total self-consisted confining po- 
tential, we can compute the Bloch wave functions and 
wave vectors by solving the eigenequation H19|) for the 
whole range of energies of interest (note that for these 
calculations the energy has to be chosen on the real axis). 
Knowledge of the wave vectors for different states allows 
us to recover the subband structure, i.e. to calculate an 
overage position xj^ of the wave functions for different 
modes a^. 



We calculate the conductance of the wire G°" = I'^ /V 
on the basis of the linear-response Landauer formula, 

where summation is performed over all propagating 
modes a for the spin a, with i^f^^t being the propaga- 
tion threshold for a-th mode. In order to visualize the 
current density we can re-write Eq. H31|) for the total 
current in the form = a J2n •^a ("-) ; where the current 
density for the mode a reads 

with and j^{n, E) being respectively the group veloc- 
ity and quantum-mechanical particle current density for 
the state a at the energy E (see Eqs. (|^ . (|^ '). and V 
being the applied voltage. 

IV. RESULTS AND DISCUSSION 

A. Hartree approximation 

To outline the role of exchange and correlation inter- 
actions we first study the magnetotransport in a quan- 
tum wire within the Hartree approximation (i.e., when 
V^^{y) is not included in the effective potential lO, and 
the spin polarization is driven by Zeeman splitting of 
the energy levels). In our calculations we use the pa- 
rameters of a quantum wire indicated in Fig. ^ and 
the temperature T = IK. With these parameters the 
effective width of the wire is ~ 400nm, and the sheet 
electron density n « 1.5 • lO^'^m"^. Figure 0] (a) shows 
the one-dimensional (ID) electron density n'^j^ for the 
spin-up and spin-down electrons in the quantum wire. 
The pronounced feature of this dependence is a 1/B- 
periodic, loop-like pattern of the density spin polariza- 
tion P„ = as illustrated in Fig. |i|(b). 

Figure 01 (c) shows the number of spin-resolved sub- 
bands as a function of B. (As the calculations are done 
for the finite temperature T, for a given magnetic field. 




B(T) 

FIG. 4; (Color online), (a) One-dimensional charge den- 
sity for the spin-up and spin-down electrons, n\j~,,n[jj in 
the Hartree approximation, and (b) spin polarization of the 

charge density, P„ = — , as a function of magnetic field 

B. (c) Total number of subbands, conductance of the spin-up 
and spin-down electrons, the total conductance G — + , 
and the filling factor in the center of the wire u{0) (note 

Gch = ^HO)- (d) The spin polarization of G and Gch- Pa- 
rameters of the wire are chosen as indicated in the caption to 
Fig. 1. Arrows in (b) indicate the magnetic fields correspond- 
ing to the magnetosubband band structure shown in Fig. |^ 
Temperature T — IK. 



we count the subbands that lie in the energy interval 
E < Ep + 4:kT, where AkT determines the energy win- 
dow beyond which the Fermi-Dirac distribution rapidly 
decays to zero). The pronounced feature of this depen- 
dence is that the number of subbands is always even, 
A'^ = 2, 4, 6, . . . , such that the spin-up and spin-down 
subbands depopulate simultaneously. The comparison of 
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FIG. 5: (Color online), (a) Spatially resolved difference in the electron densities n' [y) — n (y) as a function of B calculated 
within Hartree approximation, (b)-(e) The subband structure for the magnetic fields indicated in (a) (see also Fig. Bl(b)). 
Upper panel: The filling factor viy) for spin-up and spin-down electrons. Middle panel: the current density distribution for 
spin-up and spin-down electrons calculated according to Eq. I32II . Lower panel: magnetosubband structure for spin-up and 
spin-down electrons (solid and dashed lines correspondingly). Fat solid and dashed lines indicate the total confining potential, 
Eq. mm . for respectively spin-up and spin-down electrons. Temperature T — IK. 



Figs. H(a)-(c) demonstrates that the spin polarization is 
directly related to the magnetosubband structure: The 
polarization drops almost to zero at the magnetic fields 
when the subbands depopulate. In order to understand 
the origin of the spin polarization let us analyze the evo- 
lution of the subband structure as the applied magnetic 
field varies. Let us concentrate at the polarization loops 
in the field interval 1.3T ^ B < 2.6T when the number of 
the spin-resolved subbands = 4 and the filling factor 
in the middle of the wire 2 < 1^(0) < 4. 

Figure [3 (a) shows the spatially resolved difference in 
the electron densities n'^ (y) — n^{y) as a function of B. 
When the subband number > 4 (for B < 2.6T), the 
electron density is mostly polarized in the inner region 
of the wire. We thus concentrate first on the forma- 
tion of the compressible and incompressible strips in the 
inner region due to the upper subband. Figure |5l (b) 
shows the filling factor v''{y), current densities J^in), 
and the magnetosubband structure for the magnetic field 
B = 1.35T. This field corresponds to the case when the 
5th and 6th spin-resolved subbands just became depop- 
ulated, i.e. their bottoms are situated at > Ep + AkT 
. The 3rd and 4th subbands are separated from the 5th 
and 6th by the distance fiuJc (with lOc being the cyclotron 
frequency, HlOc ^ kT) , see Fig. |5l (b) . They are there- 
fore situated below the Fermy energy and are fully pop- 
ulated. As the results, the electron density is constant, 
which corresponds to the formation of the incompress- 
ible strip. Because of the spin-up and down subbands 
are fully filled, the corresponding electron densities are 
equal and the spin polarization of the electron density is 
zero. 

When magnetic field is raised the subbands are pushed 
up in the energy and the two highest spin-resolved sub- 



bands, following Chklovslii et al. scenario^, become 
pinned at the Fermi energy. The subband bottoms flatten 
which signals the formation of the compressible strip in 
the middle of the wire, see Fig. 111(c). When the subband 
bottoms reach the energy E k, Ep — 4fcr , the subbands 
become partially occupied. Partial subband occupation 
combined with their energy separation due to Zeeman 
interactions results in the different population for spin- 
up and down electrons. With increase of the magnetic 
field the filling factor decreases, but spin polarization in- 
creases until the subband bottoms approach ~ Ep, Fig. [S] 
(d) . This magnetic field corresponds to the maximal spin 
polarization P„ ~ 3%. With further increase of the mag- 
netic field, the subbands bottoms are pushed up above 
Ep^ which causes further decrease of the filling factor and 
diminishing screening efficiency. As the result, the width 
of the compressible strip decreases until the upper sub- 
bands become completely depopulated and the incom- 
pressible strip forms again in the middle of the wire, see 
Fig. ini(e). This is accompanied by a gradual decrease of 
the density polarization P„ to zero. The shrinkage of the 
compressible strip in the middle of the wire can be also 
clearly traced in the evolution of the current density dis- 
tribution, shown in the middle panels of Figs. (b)-(e). 
It is interesting to note that the compressible regions are 
not formed for the outermost edge states corresponding 
to the lowest subbands A^ = 1 and 2. This is because 
that in the field interval under study the extension of the 
wave function is larger than the width of the compress- 
ible strip predicted by the Chklovskii et al. theory^. The 
onset of the formation of the compressible strips can be 
seen in Fig. (e) for B = 2.5T. Note that the effect of 
the formation/non- formation of the compressible strips 
in quantum wires was discussed in details by Suzuki and 
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y (nm) 

FIG. 6: (Color online) Lower panel: A closeup of the spin- 
up and spin-down subbands A'^ = 3, 4 (solid and dashed lines 
respectively) for the magnetic field B = 1.5T, when the po- 
larization of the conductance Pq is negative. Upper panels; 
current density distribution for the spin-up and spin-down 
electrons at the energy E indicated in the lower panel. At 
this energy, there is one propagating state for spin-up elec- 
trons and three propagating states for spin-down electrons. 
Left inset shows the derivative of the Fermi-Dirac distribution 
determining the weight of the contributions from the current- 
carrying states to the total current density at the given energy, 
see Eq. CT . 

Ando for the case of spinless electronsi^. 

The described above picture of evolution of the density 
polarization qualitatively holds for all other polarization 
loops. We stress that in all the loops only two upper, par- 
tially occupied spin-resolved subbands contribute to the 
spin polarization, whereas remaining subbands are fully 
(and equally) populated and thus do no contribute to 
the total spin polarization. When magnetic field exceeds 
B = 2.6T, only two subbands survive in the quantum 
wire. With further increase of magnetic field the upper 
(spin-up) subband gradually depopulates and the density 
polarization P„ grows linearly until it reaches 100% when 
only the spin-down subband remains in the wire. 

It should be also stressed that within Hartree approxi- 
mation two outermost spin-up and spin-down edge states 
are not spatially polarized (i.e. they are situated at a 
practically same distance from the wire edges, see Fig. 

13. 

Figures 01 (c),(d) show the conductance G"^ for spin-up 
and spin-down states and its relative spin polarization 

T+Gi ■ "^^^ ^Pi'^ polarization Pq follows a similar 
behavior as the density polarization P„ with one subtle 
difference. Namely, the density polarization P„ is always 
positive because spin-up states always lie in energy be- 



low the corresponding spin-down states, and, therefore 
{y)—n^{y) > 0. In contrast, the spin polarization of the 
current, after reaching zero, does not immediately raises 
as the magnetic field increases, but, instead, becomes 
negative before raising again. Note that this is accom- 
panied by a small (but noticeable) increase of the total 
current (at B ~ 1.5T, ST, see Figs. 111(c)). This effect 
can be traced back to the self-consistent band structure 
as explained below. FigurejBlshows a closeup of the upper 
subbands = 3,4 for the magnetic field B = 1.5T, i.e. 
when the current polarization is negative. Because the 
spin-up/down subbands are not flat, for certain energies 
E < Ep the upper (spin-down) subband can give rise 
to several propagating states, whereas the lower (spin- 
up) subband corresponds to only one propagating state, 
see Fig. (HI According to the Landaulet formula ^Al\ all 
propagating states contribute equally to the total cur- 
rent. Because of this and due to the fact that the spin- 
down subband is situated closer to the Fermi energy, the 
total current for the spin-down electrons is larger than 
the current for the spin-up ones. This explains the neg- 
ative spin polarization of the current and the increase of 
the total current at the magnetic fields just above the 
subband depopulation. We are not aware of the discus- 
sion of this effect in the current literature. The available 
experimental data, see e.g.. Fig. 2 of Ref. IssI showing a 
nonmonotonic dependence of the conductance of a quan- 
tum wire as a function of magnetic field, are consistent 
with the predicted behavior of the total current. Note 
that this feature in the conductance also survives within 
the DFT approach (see below. Fig. [TJ. 



Figures01(c) and (d) also show the conductance and its 
spin polarization calculated according to the Chklovskii 
et al. prescription^; Gch = x'^^'^)' '^ith v[{)) being the 
filling factor in the center of the wire. Gch follows the 
exact conductance rather good, but does not recover the 
steps in the conductance related to the subband depopu- 
lation (see Ref. fSi] for a related discussion). Gch does not 
also reproduce the increase of the current and the neg- 
ative conductance polarization discussed above because 
these features are related to the quantum-mechanical 
band structure. 



As we mentioned before, the Hartree approximation 
predicts that spin-up and spin-down subbands depopu- 
late simultaneously and thus the conductance drops in 
steps of 2e? jh as the magnetic field increases. This is in 
strong disagreement with the experimental observations 
that demonstrate that the subbands depopulate one by 
one such that the conductance decreases in steps of jh. 
We will show in the next section that accounting for the 
exchange and correlation interactions leads to qualita- 
tively new features in the subband structure and brings 
the theory to a close agreement with experiment. 
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FIG. 7: (Color online), (a) One-dimensional charge density 
for the spin-up and spin-down electrons, i 



nj^ calculated 

within DFT-I-LSDA, and (b) spin polarization of the charge 
density, Pn = j^, as a function of magnetic field B. 



(c) Total number of subbands, conductance of the spin-up 
and spin-down electrons, the total conductance G = + , 
and the filling factor in the center of the wire, (d) The spin 
polarization of the conductance Pa- Parameters of the wire 
are chosen as indicated in the caption to Fig. 1. Arrows in 
(b) indicate the magnetic fields corresponding to the mag- 
netosubband band structure shown in Fig. |S| Temperature 
r= IK. 



B. Density functional theory in the local spin 
density approximation 

Figures lTlSl show the electron density, conductance and 
subband structure for the quantum wire calculated us- 
ing DFT within LSDA, Eqs. g!)-®. Utilization of 
the DFT+LSDA leads to several major quantitative and 
qualitative differences in comparison to the Hartree ap- 
proximation. First, the spin polarization of the electron 



density also shows a pronounced l/_B-periodic loop-like 
pattern. However, for the given magnetic field B the spin 
polarization in the quantum wire calculated on the basis 
of DFT is of the order of magnitude higher in compari- 
son to the Hartree approximation. Second, the magne- 
tosubbands depopulate one by one, and the conductance 
decreases in steps of e'^/h (not in steps of 2e'^/h as in 
the case of Hartree approach when the spin-up and spin- 
down subbands depopulate simultaneously). Third, the 
outermost edge states become spatially polarized (sepa- 
rated), which is in the strong contrast with the Hartree 
approximation, where they are situated practically at the 
same distance from the wire boundary. 

In order to understand the effect of the exchange and 
correlation interactions on the evolution of the mag- 
netosubband structure, let us now concentrate on the 
same field interval as studied in the previous section, i.e. 
when the number of the subbands lies in the interval 
3 < iV < 4 and the filling factor in the middle of the wire 
2 < 1/(0) < 4, see Fig. [7| We start from the magnetic 
field B = 1.4T, where the spin polarization of the density 
is minimal. Similarly to the case of Hartree approxima- 
tion (Fig. El(b)), this corresponds to the case when 5th 
subbands just became depopulated as shown in Fig. |S| 
(b). However, in contrast to the Hartree approximation, 
where the Zecman interaction is not strong enough to 
cause any significant spin polarization, in the present case 
the exchange interaction leads to a non-negligible spin 
polarization near the boundaries of the wire (P„ ~ 7%). 
For this magnetic field the number of subbands is even, 
and the spin-up and spin-dow subbands are fully filled in 
the center of the wire. As the result, the electron den- 
sities are constant, which corresponds to a formation of 
the incompressible strip in the center of the wire. Be- 
cause the spin-up and spin-down subbands are equally 
occupied, spin polarization in the center of the wire is 
zero (Fig. |a(a)). 

When the magnetic field increases, all the subbands 
are pushed up in energy, and 4th subband gets pinned at 
Ep near the boundary of the wire, forming the compress- 
ible strip. Fig. IHl(c). With increase of magnetic field, the 
compressible strip extends to the center of the wire, com- 
pare Figs. IHl(b) and (c). Note that in the Hartree case the 
separation between the subbands caused by the Zeeman 
splitting is small (^ kT) and hence both the subbands 
are pinned at Ep (see Fig. |Sl(c)-(c)). In contrast, in the 
present case only one of the subbands is pinned at Ep 
because the subband separation is determined by the ex- 
change interaction whose magnitude can be comparable 
to hojc- 

Figure|Sl(d) shows the subband structure for the mag- 
netic field B = 1.9T when the spin polarization of the 
electron density is maximal. In this case 4th subband 
is about to be depopulated and all the remaining sub- 
bands (two spin- up and one spin down) lie below Ep. 
They are therefore fully populated (2n|^ « "id)' which 
corresponds to the calculated polarization P„ ~ 33%. 

When magnetic field is increased by only 0.05T, the 
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FIG. 8: (Color online), (a) Spatially resolved difference in the electron densities (y) — n^iy) as a function of B calculated 
within DFT+LSDA. (b)-(e) The subband structure for the magnetic fields indicated in (a) (see also Fig. |7|(b)). Upper panel; 
The filling factor u{y) for spin-up and spin-down electrons. Middle panel: the current density distribution for spin-up and 
spin-down electrons calculated according to Eq. Lower panel: magnetosubband structure for spin-up and spin-down 

electrons (solid and dashed lines correspondingly). Fat solid and dashed lines indicate the total confining potential, Eq. 1101 . 
for respectively spin-up and spin-down electrons. Temperature T = IK. 
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FIG. 9: (Color online). Spatial separation between the out- 
ermost spin-up and spin-down edge states as a function of 
magnetic field B. (The separation between the edge states 
is extracted from the corresponding current distribution, see 
Fig. |H| middle panel) . The number of subbands and the elec- 
tron density spin polarization P„ from Fig. Q is shown for 
comparison. 



density spin polarization drops by ~ 10%, and the sub- 
band structure experiences dramatic changes, see Fig. |S| 
(c). In particularly, the spatial separation between the 
outermost spin-up and spin-down states collapses from 
^ 20 nm to zero, as shown in Fig. |51 The explanation 
of this remarkable effect is based on the fact that the 
electrostatic energy is dominant for the system at hand^. 
This is illustrated in Fig. llOl which compares the electron 
densities and the magnetosubband structure in a quan- 
tum wire calculated within the Hartree and DFT approx- 
imations for some representative value of the magnetic 
field. As expected, the total electron density is practi- 
cally the same in the both approximations. At the same 



time, the magnctosubbands and the spin-up and spin- 
down densities vary significantly between them. It is also 
interesting to note that the magnetic fields correspond- 
ing to the depopulation of even subbands TV = 2, 4, 6 . . . 
are practically the same with and without accounting for 
the exchange and interaction terms, compare Figs. 0] (c) 
and[7|(c). The dramatic changes in the subband struc- 
ture at i? ~ 1.95 T can be explained as follows. At 
B ^ 1.9T the electron density near the edge of the wire 
is dominated by spin- up electrons, see[7|(d), the upper 
and middle panels. When magnetic field is raised, 4th 
subband practically depopulates, and 3rd (spin-up) sub- 
bands is pushed up in energy. As the result, the density 
of the spin-up electrons associated with this subband is 
redistributed towards the center of the wire. However, 
this small change in the magnetic field can not affect 
the total density. Because of this, the density of the re- 
maining electrons has to be adjusted to keep the total 
density unchanged. This can be done only if the spin- 
down electrons associated with the subband N ^ 2 are 
redistributed towards the edge of the wire. As a conse- 
quence of this redistribution, the densities of the spin-up 
(1st subband) and spin-down (2nd subband) electrons 
near the wire edge become approximately equal and so 
does the total confining potential V'^{n), Eq. 1(1171) . The 
latter results in the absence of the spatial separation for 
the outermost edge states = 1 and N = 2. Note that 
the effect of a collapse of the spatial separation between 
the outermost edge states is related to the features of the 
quantum-mechanical band structure, and hence is absent 
in the Thomas-Fermi approximation^. This effect can 
be utilized in spintronics devices operating in the edge 
state regime for injection of different spin speciesi^. 

The outermost spin-up and spin-down edge states re- 
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FIG. 10: (Color online). The filling factor, current densities and the magnetosubband structure (upper, middle and lower 
panels correspondingly) calculated within Hartree and DFT approximations for two quantum wires with different distances 
between the gates, (a) a = 500 nm, and (b) a = \\im. Remaining parameters of the wire are chosen as indicated in the caption 
to Fig. 1. Solid and dashed lines correspond to the spin-up and spin-down states. Fat solid and dashed lines indicate the total 
confining potential, Eq. dU, for respectively spin-up and spin-down electrons. Temperature T = IK. 



mains spatially degenerate up to the magnetic field B ^ 
2.25 T, see Fig|Hl(a) and Fig. El The spin polarization of 
the electron density P„ gradually decreases in the range 
1.9T< B < 2.8T. This decrease is related to the gradual 
depopulation of the 3rd (spin-up) subband. At B ^ 2.8T 
this subband practically depopulates, P„ reaches its min- 
imum, and the incompressible strip is again formed in the 
middle of the wire (Fig. |Hl(f))- With further increase of 
the magnetic field, 2nd (spin-up) subband gets pinned 
to Ep, and P„ gradually increases imtil it reaches 100% 
at the magnetic field when the spin-up subband depopu- 
lates. 

As in the case of Hartree approximation, the evolution 
of the magnetosubband structure within DFT described 
above qualitatively holds for all other polarization loops. 

Figure [3c) shows the conductance for spin-up and 
spin-down electrons , , the total conductance G = 
G^ -|- G''-, the filling factor in the middle of the wire z^(0), 
and the spin polarization of the conductance. The total 
conductance G{B) decreases in steps of e^//i closely fol- 
lowing the depopulation of the magnetosubbands as B 
increases. Note that the magnitude of G{B) in plateau 
regions when N > 2 shows slight increase in compar- 
ison to the corresponding value of Ne'^/h. This effect 
has the same origin as in the case of Hartree approxi- 
mation (see Fig. Eland a related discussion in the text). 
For N < 2 this effect becomes much more pronounced 
in comparison to the Hartree approximation. This is be- 
cause for magnetic fields corresponding to < 2, the 
separation between bottoms of spin-up and spin-down 
subbands due to the exchange interaction exceeds 8kT. 
Because the subbands are not flat, the spin-down sub- 
band (which is pinned to Ep) gives rise to several states 
propagating in the bulk of the wire as discussed in the 
previous section, whereas the spin-up subband (whose 



bottom lies well below Ep) corresponds to only one prop- 
agating state situated near the wire edge. 

Note that the propagating states giving rise to the con- 
ductance for N < 2 are the Bloch states of an infinite 
quantum wire. In a typical experimental condition, a 
long quantum wire is connected to a much wider region 
of 2DEG^. The edge states in the region of 2DEG are 
coupled only to the edge states in the wire. As the results, 
the measured conductance for N < 2 does not exhibits 
the increase over the plateau values of Ne^ /}&. 

Finally, we note that our analysis of the spin polar- 
ization and evolution of magnetosubbands in quantum 
wires was concentrated on a representative wire with the 
distance between the gates a = 500 nm and the sheet 
electron density n w lO^^m"^. We would like to empha- 
size that all the results presented here qualitatively hold 
for wires of arbitrary widths and electron densities. This 
is illustrated in Fig. 1 101 for the case of two quantum wires 
with different distances between the gates, a = 500 nm 
and a = 1/im, which shows practically identical subband 
structure as well as electron and current densities distri- 
butions. 



V. CONCLUSION 

In the present paper we provide a quantitative de- 
scription of the structure of edge states in split-gate 
quantum wires in the integer quantum Hall regime. 
We start with a geometrical layout of the wire and 
calculate self-consistently quantum-mechanical magneto- 
subband structure and spin-resolved edge states where 
electron- and spin interactions are included within the 
density functional theory in the local spin density ap- 
proximation (DFT+LSDA). 
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We develop an effective and stable numerical method 
based on the Green's function technique capable of deal- 
ing with a quantum wire of arbitrary width in high per- 
pendicular magnetic field. The advantage of this tech- 
nique is that it can be directly incorporated into magne- 
totransport calculations, because it provides the eigen- 
states and wave vectors at a given energy, not at a given 
wavevector (as conventional methods do). Another ad- 
vantage of this technique is that it calculates the Green's 
function of the wire, which can be subsequently used in 
the recursive Green's function technique widely utilized 
for magnetotransport calculations in lateral structures. 

We use the developed method to calculate the self- 
consistent subband structure and propagating states in 
the quantum wires in perpendicular magnetic field. We 
discuss how the spin-resolved subband structure, the cur- 
rent densities, the confining potentials, as well as the 
spin polarization of the electron and current densities 
evolve when an applied magnetic field varies. We demon- 
strate that the exchange and correlation interactions dra- 
matically affect the magnetosubbands in quantum wires 
bringing about qualitatively new features in comparison 
to a widely used model of spinless electrons in Hartree 
approximation. These features can be summarized as 
follows. 

(a) The spin polarization of the electron density shows 
a pronounced l/i?-periodic loop- like pattern, whose pe- 
riodicity is related to the subband depopulation. For a 
given magnetic field B the spin polarization in the quan- 
tum wire calculated on the basis of DFT+LSDA is of the 
order of magnitude higher in comparison to the Hartree 
approximation (where the spin polarization is driven by 
the Zeeman interaction only). 

(b) The magnetosubbands depopulate one by one, and 
the conductance decreases in steps of e^/h (not in steps 
of 2e'^ /h as in the case of Hartree approach when the 
spin-up and spin-down subbands depopulate practically 
simultaneously) . 

(c) The outermost spin-up and spin-down edge states 
become spatially polarized (separated), which is in the 
strong contrast to the Hartree approximation, where they 
are situated practically at the same distance from the 
wire boundary. We also find that the spatial separa- 
tion between the outermost edge states disappears in the 
range of magnetic close to filling factor 1/ ~ 3 and then 
is restored again when the magnetic field is raised. This 
effect can be utilized in the spintronics devices operat- 
ing in the edge state regime for injection of different spin 
speciesi^. 

Recently, the structure of edge states around quantum 
antidots has been the subject of a lively discussion^S. 
Even though the method developed in the present paper 
applies to quantum wires, it is reasonable to expect that 
for sufficiently large antidots (when the single particle 
level spacing A is smaller than kT) the present approach 
can also provide information on the edge state structure 
around the antidots. 

A direct probe of spin polarization of electrons in quan- 



tum dot edge channels using polarized photolumines- 
cence spectra has been recently reported by Nomura and 
Aoyagi^. Their method opens up a possibility for a di- 
rect probing of the electron density spin polarization in 
quantum wires, such that the results presented in our 
study (in particularly the spin polarization shown in Fig. 
Ob) and Fig. |HIa)), can be directly verified in the exper- 
iment. 
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APPENDIX A: EXCHANGE AND 
CORRELATION POTENTIALS IN THE LOCAL 
SPIN DENSITY APPROXIMATION 

In this Appendix we provide explicit expressions for 
the exchange and correlation potentials entering the DFT 
effective potential The exchange and correlation en- 
ergies for 2DEG used in Eq. |(HJ) are given by Tanatar 
and Ceperley (TC)S1. The exchange energy reads 



-Ry" 



4V2 
37rr» 



(i + e)^ + (i-e)^ 



(Al) 



where is the dimensionlcss density parameter which 
is defined in terms of the effective Bohr radius (ap- 
propriate for a material with the effective electron mass 
m* = rrieffrne, and the dielectric constant e = SrSo), 



£r 



meffniee 

where Bohr radius ao = 0.529 • 10^*^ 

Ry* 



-ao, 



(A2) 

m. The factor 



'-^Ry {IRy = 2.17989 • IQ-^^J) 



generalizes TC results for the case of an arbitrary effec- 
tive electron mass m* and relative dielectric constant £r , 
and converts TC expressionsSi into SI units. The corre- 
lation energy for the unnpolarized case (^ — 0) and for 
the fully polarized case (^ = 1) is approximated in the 
formal. 



EcoriO = -Ry*Co 



1 + Ciw 



1 + Ciw + C2uP + C3W^ ■ 



(A3) 



where w = and the coefficients Co,Ci,C2,C3 are 

tabulated below For the case of an intermediate polar- 
ization, < ^ < 1, the correlation energy can be inter- 
polated between the nonpolarized and the fully polarized 
cases following the receipt of von Barth and Hendin2^i2i 

EcoriO - EcoriO) + fiO [Ecoril) - EcoriO)] (A4) 



with fiO = 



(i + 0i + (i-0i-2 

23/2 _ 2 
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unnpolarized case 
(C = 0) 


fully polarized case 
(C = l) 


Co 


-0.3568 


-0.0515 


Ci 


1.13 


340.5813 


C2 


0.9052 


75.2293 


C3 


0.4165 


37.0170 



TABLE I: Tabulated coefRcients Co,Ci,C2,C:i (see text for 
detail). 



V^oriO = Vcor{0) + fiO [Vcoril] " V,o.(0)] 

^cor(e = oore = i) = 



(A6) 



-RyCo 



1 -I- diw -t- d2w'^ + dsw^ + diw^ 



{1 + Ciw + C2W^ + Csw^r 



Taking the functional derivatives © using the above 
expressions for the exchange and correlation energies 
(|Aip . l|A3|l we arrive to the following expression for the 
exchange potential Vex^, Vexi, and for the correlation po- 
tentials used in Eq. 



^/2 



ex] 



Ve 



4 eqe 

P 
V2 



^VP [(1 + 0^ + (1-0^ 



ex I 



VpI [(1 + 0* + (1-0^ 



4 eo^rTr^/^ 



2PT 



(A5) 



where /(O is given by Eq 



^C2 + Cf), d3 = (1^3 + 1^1 C2), d 



41, and di 

4 : 



2Ci, d2 - 
IC1C3. 
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